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Abstract 

7^ , 1 The anomalously large like-sign dimuon charge asymmetry in semileptonic 5-hadron decays recently 
measured by the DO Collaboration may be hinting at the presence of CP-violating new physics in the 
mixing of Bg mesons. It has been suggested that the effect of a nonstandard spin-1 particle lighter than 
CD the b quark with flavor-changing couplings to h and s quarks can reproduce the DO result within its 



00 



one-sigma range. Here we explore the possibility that the new particle also couples to charged leptons 
i = e, n and thus contributes to rare b ^ s processes involving the leptons. We consider in particular 
fSj ' constraints on its couplings from existing experimental data on the inclusive B — )• Xs£~^i~ and exclusive 
B — )• K^*^i^i^ decays, as well as the anomalous magnetic moments of the leptons. We find that there 
' is parameter space of the particle that is allowed by the current data. Future measurements of these B 
^ I transitions and rare decays of the Bg meson, such as Bg — s- {4>,'r],r]')£'^i^ and Bs — s- at LHCb and 

i next-generation B factories can probe its presence or couplings more stringently. 
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I. INTRODUCTION 



The DO Collaboration has recently reported a new measurement of the like-sign dimuon charge 
asymmetry in semileptonic 6-hadron decays jl| which disagrees with the standard model (SM) 
prediction by about three standard deviations. Although this finding still needs to be confirmed 
by future experiments, it might be a clue to the presence of unexpectedly substantial CP- violating 
new physics in the mixing of Bg mesons. Subsequently, we have suggested as one of the possible 
scenarios for the new physics that the DO result could be attributed to a nonstandard spin-1 particle 
lighter than the b quark, with fiavor- changing couplings to the b and s quarks 0]. Specifically, we 
showed that the effect of the new particle can lead to a prediction which is consistent with the new 
data within its one-sigma range. 

New-physics scenarios involving nonstandard spin-1 particles with masses of a few GeV or less 
have been discussed to some extent in various other contexts in the literature. Their existence 
is generally still compatible with currently available data and also desirable, as they may offer 
possible explanations for some of the recent experimental anomalies and unexpected observations. 
For instance, a spin-1 boson having mass of a few GeV and couplings to both quarks and leptons 
has been proposed to explain the measured value of the muon g—2 and the NuTeV anomaly 
simultaneously |i3i]. As another example, 0{MeV) spin-1 bosons which can interact with dark matter 
as well as leptons may be responsible for the observed 511-keV emission from the galactic bulge 0]. 
If its mass is at the GeV level, such a particle may be associated with the unexpected excess of 
positrons recently observed in cosmic rays, possibly caused by dark-matter annihilation In 
the context of hyperon decay, a spin-1 boson with mass around 0.2 GeV, flavor-changing couplings 
to quarks, and a dominant decay mode into can explain the three anomalous events of 

S+ — )■ pfi'^fi~ reported by the HyperCP experiment several years ago jsj-lsj . Although in these few 
examples the spin-1 particles tend to have suppressed couplings to SM particles, it is possible to 
test their presence in future high-precision experiments |J-l9|. It is therefore also of interest to see 
if the light spin-1 boson that could be responsible for the anomalous DO result contributes to some 
other 6-meson processes, perhaps with detectable effects, which we will attempt to do in this paper. 

Here we consider the possibility that the new spin-1 particle, which we shall refer to as X, has 
flavor- conserving couplings to the electron and muon, besides its flavor- changing couplings to the 
b and s quarks. Accordingly, it can contribute to a number of rare 6 — )■ s transitions involving the 
charged leptons, via the quark-level process b — j- sX* — )■ si~^i~, where i = e,fi. In particular, we 
will deal with the impact of X on the inclusive 6-meson decay B — )■ XJ~^£~ and the exclusive ones 



B —7- K'^*^£'^i~ and B^ — )■ (pi^i" , all of which have been observed 10l-ll4|. as well as the leptonic 
decay B^ — > . Using the existing experimental information on the 6- meson decays as well as 
the recent data from DO on Bg mixing, we will explore constraints on the couplings of X. Since 
X has flavor-conserving interactions with the leptons and hence contributes to their anomalous 
magnetic moments, we will also take their measured values into account. 

In the following section, we write down the relevant Lagrangians and derive the amplitudes for 
the processes of interest. In Sec. IIIIl we present our numerical results for the constraints on the X 
couplings and provide some predictions which may be tested at LHCb or future B factories. We 
conclude in Sec. |IV] and collect some of the formulas in appendixes. 
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II. INTERACTIONS AND AMPLITUDES 



We adopt a model-independent approach, assuming in addition that X carries no color or electric 
charge and its couplings to the fermions have both vector and axial-vector parts. The Lagrangian 
for its flavor- changing interactions with the h and s quarks is then ^ 



bsX 



-si^,{gvs- 9Asl^)bX^ + H.c. 



(1) 



where gy^ and g^^^ parametrize the vector and axial-vector couplings, respectively, and in general 
can be complex, which would be new sources of CP violation. Moreover, the flavor-conserving 
couplings of X with a charged lepton ^ is described by 



tx 



(2) 



where gyg^ and (7^^ are real parameters because of the hermiticity of C^x- We will study only 
the ^ = e and cases due to lack of the relevant data for £ = r at present. In the absence of 
model specifics, gy^^Ae gy^j^^Afi necessarily related. In principle, X can have additional 

interactions, flavor-conserving and/or flavor- violating, with other fermions which are parametrized 
by more coupling constants. We assume that these additional parameters already satisfy other 
experimental constraints to which they are subject, but which we do not cover in this paper. 

Together C^^^x '^ex generate the contributions of X to the above-mentioned 6 — )■ s transi- 
tions. For the inclusive decay B — )■ Xs£~^£~, the resulting amplitude is 



M 



X 

b^sU 



2g^^m^s[{m^ - m^)gy^ + (m^ + m^)gj^^-i^]hh^£ 



m 



X 



(g2 



m 



X 



+ i^x^x) 



(3) 



where q = +Pe- is the combined momentum of the dilepton and Tj^ the total width of X. We 
remark that the presence of the dependence in the denominators distinguishes this new-physics 
scenario from others involving heavy particles, which would induce four-fermion operators with 



coefficients independent of g^, as have been studied in the literature [15 
is well known and given by 



0— >-s 



16| . The SM counterpart 



M 



SM 



V2 



— ^ q^ -sa^^ (m, P^ + m^Pj & h/ 



(4) 



where a = e^/(47r) and Gp denote the usual fine-structure and Fermi constants, respectively, 
Vki are Cabibbo-Kobayashi-Maskawa matrix elements, Cy^ are Wilson coefficients evaluated at 
a scale /i ~ m^, and -Pl r = |(1 -FTs)- From the sum of the SM and X-induced amplitudes follows 
the decay rate 



X 

b^sU ' 



(5) 
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where 



^SM _ 
b^sU 



768 vr^ ml 



\Cf\' + \Cf,\'){ml + 2q^) + 



eff I 



mt 



{2ml + q') 



+ 12Re{Cf*Cf) 



[mb-q ) 



(6) 



^^'^^ 768 vrS ml 



dq^ { Re 



^ /^cff* ^ ^cff* ^ \ S'ys + 9 As 



A 



[ml + 2g2) 



+ 6Re( nCf* ^^^^ )^7v,m^ " 



X 



with and having been set to zero and 

2v^7r 



A, 



ts^ tb 1 



K 



A 



X 



2 2 I -n 

q — mj^ + tL x^x • 



(7) 



(8) 



Evidently T^^^^ contains not only the X-induced amplitude alone, but also its interference with 



the SM one. 

For B — K£~^£~, the amplitude follows from the effective Hamiltonians "H^^^^ which yield 
Eqs. ([3]) and (jll). Adding the X and SM contributions, we can write it as 



—aGp \ 
2V2n 

where q = + = Pb - Pk, 



A [pb + PKYh,^ + [C [pb +PkY + D q']h,l,^ 



(9) 



A 



D 



-cff , '^9vs9vt\p^ , 2m^CfFT 



Ax 



+ 



m^ + mj^ 



C 



^off _ ^ 9vs 9ai j 



A 



X 



C- 



•'K 
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m 



B 



m\ K Qy^ g^^ [F^ m\ + Fq (g^ _ m\)] 



m\ q^ 



A 



(10) 



X 



with -Fq 1 T denoting the form factors for the B ^ K matrix elements of the h ^ s quark operators 
in T-C^'^i^ and being defined in Appendix [Aj One can see that the B — )■ Ki^i~ amplitude is 
independent of g^^g. The corresponding decay rate is given in Appendix iBl 

It is worth mentioning here that the i^^-meson decay B^ — t- P£^£^ involving a pseudoscalar 
meson P containing an ss component in its quark content, such as 1] and rj', also has an ampli- 
tude independent of gj^^. Consequently, the observation of B^ — )■ P£'^£^ will provide additional 
information on the X contributions involving the products gy^ gy^ and gy^ g^^. 

For B — 7- K*£^£~, the amplitude from the SM and X contributions can be expressed as 



M{B^ K*£+£-) 



2V2 



TT 



A e^,^^ -tCe; + tVe*-q {p^ + Pk^)^ 



+ 



Se^^^^e*''p-BPK^-iFel + ige*-q{pB+PK^), + tne*-qq^Iri,(] , (H 



where q = + j9^„ = Pb - Pk* > 



A = cr + — — + 



^ Ax Jm^ + nij^, g2 

\ 2 2 



C = {cf + ^^f^y,{m, + m^,)+2m,CfT, 



2 ^2 



-yeff I ^ 9 As 9vi \ ^2 , o_ ^off ( ^2 , ^3 



9 "I \ r ^'iiL^Ky-j o -| 5 ^ 



K* 



F - / - *1^Xi1m\ F - /^r'^ff - ^ ^-^^ ^^f- \A{m \ 



cff S'yls fi'Af \ ^2 



the form factors V ^ ^012; 2 3 ^^e B — )■ A'* transition being defined in Appendix |Al The 

corresponding squared amphtude and decay rate are given in Appendix iBl 

From the squared amphtude for B — t- K*t^l~ ^ one can arrive at two additional observables 
which have been measured besides the branching ratio. They are the longitudinal polarization 



fraction and the lepton forward-backward asymmetry Apg, which are defined from 18|, |l9 



1 d^T(B ^ K*i+i~) 3, , , 3, , , 

^ ^ - .(l-cos2^)F^ + -(l + cos2^)(l-F^) 



dT(B ^ K*£+i-)/dq^ dq'^d{cos9) 4 

+ AfbCos^, (13) 

where 6 is the angle between the directions of B and i~ in the dilepton rest frame. Since these 
observables are ratios of squared amplitudes, their dependence on the hadronic form-factors is partly 
canceled, which reduces the theoretical uncertainties associated with the form factors. Espec ially 



y4pB is predicted with good precision in the SM to have a zero-crossing point at ~ 4 GeV^ 18 
which makes this asymmetry very sensitive to the signals of new physics that can shift the location 
of the point. In Appendix [B] we have written down and Ap^ in terms of A,C, . . . ,1-1. 

The amplitude for the contribution of X to -B^ — i- i^i" can be obtained from Eq. (|3]) after 
making the approximation q^ = m^^ in the center-of-mass frame of 6s and neglecting the part. 
Thus, using the matrix elements (0|s7'^6|5^) = (0|s6|i?^) = 0, (0 1 57^756! 5^ (p)) = —ifssP^^ and 
{0\s'y^b\B^) = ifssmlj (m^ + m^), we arrive at (sj 



Alf^,,,_ = -^-llM.^Al3Mj!hl^^i, (14) 



mx 



Since {0\sa^'''b\B^) = {o\sa^"'-f^b\B^) = 0, the SM yields 



<SM iaG-p\ /s^m^ ^ - 



= ^i^^C-W- (15) 
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The sum of the two results in the decay rate 



r(4^£+r) 



167r3 



rics , ^ 9 As 9ae 

•-"10 ^ 



m 



X 



m\ — Amj . (16) 



Hence rl does not probe either gy^ or gy^^. 

The X interaction with the lepton ^ as described by C^^x Eq. (|2]) affects the anomalous 
magnetic moment of i at one loop. The X contribution to can be expressed as [g], Eo| 

2 

ofK) = ^^(^?w/y(0+^7A./AW) , (17) 



where r = m\/m\, 



^ — „ , , , —Ax + — (1 + 2r)x^ 



/^(r) = / dx- — ^, /^(r) = / dx " 2 ■ (18) 

Q 1 — X + rx^ Jo I — X + rx^ 



Since the anomalous magnetic moments of the electron and muon, and a^, have been measured 
precisely, we need to take their constraints into account. 



III. NUMERICAL ANALYSIS 
A. Constraints 

To obtain the first set of constraints on the contributions of X, we employ the data on the 
inclusive decay B — )■ Xsi~^i~ . The BaBar and Belle Collaborations have measured its branching 
ratios B for different ranges of the squared dilepton invariant mass, q^. We take their results for the 
low- and high-g^ ranges 1 GeV^ < Q'^ < 6 GeV^ and > 14.4 GeV^, respectively. The numbers 



from BaBar, B^°^ = (1.8 ± 0.7 ± 0.5) x 10"^ and B^^f = (5.0 ± 2.5l[]:f) x 10~^ [10|, and from 
Belle, B^^ = (1.49 ± 0.501°;^^) x 10"^ and B]^^f = (4.2 ± 1.2+°;^) x 10"^ [11], average out to 

-^5(5 ^ XJ+r) = (1.6 ± 0.5) X 10-6 , B^lf{B XJ+r) = (4.4 ± 1.2) x 10"^ . (19) 



We also need the SM predictions |21[ | 

13sm{B X,e+e-) = (1.64 ± 0.11) x lO"*' , B^i^{B -> = (1.59 ± 0.11) x 10"^ , 

Bl^\B ^ X,e+e-) = 2.09 x 10-^(1^°;^) , BS{B -> = 2.40 x 10-^(ll°:i) . (20) 

Upon comparing these experimental and SM values, we can require that the X contributions, from 
r^^^-^ in Eq. ([7]), satisfy 

- 5 X 10-^ < i3^"(fi ^ XJ+r) < 4 X 10"^ , < Bx^\B X^tr) < 3.5 X 10"^ , (21) 

where 

BxiB^XJ^r) = T^T^_,^,, (22) 
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with being the B hfetime. It is clear from Eq. ([7]) that Bx{B — )■ Xgi'^i ) contains both the 
X-mediated amplitude and its interference with SM one. 



Numerically, we use = +Tbo) ~ 1-582 ps, the average of the B+ and B° lifetimes m 



a = tt(mj = 1/133, Gp = 1.166 x 10"^ GeV-^, = V;Xb = -0.0405 + 0.0007? j23|, and 
the lepton and meson masses from Ref. [2^, as well as the Wilson coefficients C^^ = —0.304, 



Cg = 4.211 + Y{q ), and C^q = —4.103 from Ref. [16|, the expression for the complex function 



y(g^) given therein. These input parameter values are also used in the rest of the paper. Since 
'^b^'^e -^l- © behave as ml if gets large, they have sizable uncertainties depending on the 
choice of value. In our numerical treatment of T^^^^, we set = 4.5 GeV, as its application in 
i3gj^ = T^Tl^^^^ leads to B]Z{B X,^+r) = 1.70 x'lO-^ and Bs'i'{B XJ+r) = 1.66 x 10"^ 
which are compatible with the predicted ranges in Eq. (1201) from more refined calculations. 

The second set of constraints comes from the data on B ^ K^*H+i-. The BaBar, Belle, 



and CDF Collaborations have measured several different observables in these decays |12|-|14l]. We 
will choose the branching-ratio results provided by Belle and CDF, as they are available for the 
specific ranges for which the most recent predictions in the SM with detailed estimates of the 
uncertainties are also available. Their numbers are listed in Tables |T] and [TTl respectively. In view 
of the currently sizable experimental and theoretical errors, we will ignore the numerical effects of 
B^-B^ and e-yU differences on these processes. Furthermore, we will take the more precise of each 
pair of experimental values in Table [T] and then compare it to the corresponding SM number in 
Table [III in order to estimate the allowed range of the X contribution Bx ■ Accordingly, we can 
impose the limits 

- 0.7 X 10-^ < Bx{B^ ^^^^i,2g[i,6]GeV^ ^ 0-4 X 10-' , (23) 
-3x10-^ < i3^(5->ir*£+r)^,^j^g,^^^, < 0.5xl0-^ 
-0.5x10-' < Bx{B^K*tr)^,^^^^^^^^^^^^^, < 0.7x10-', (24) 



-0.1x10-' < Bx{B^K*tr)^,^^^^^^, < 1.1x10 



-7 
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where 



Bx{B^K^*k^r) = tbTx{b ^ i&k+r) , (25) 



with T x{B — 7> K^*H~^i ) being the rates of -B — ?► K^*H~^i , whose formulas are given in Ap- 
pendix [Bl minus the purely SM part. 

The K* longitudinal polarization fraction and lepton forward-backward asymmetry A-p ^ in 
B — i- K*i~^i~ have also been measured by BaBar, Belle, and CDF. Although the BaBar 



12 



and Belle 13| data on A^^ (^B — )■ K*i'^i ) exceed the SM expectation, the most recent CDF 
measurement IJ] is consistent with it. It is expected that more definitive information on this and 



other observables will be available with the upcoming results from LHCb in the near future. Since 



most of the current data on and Ap-^ still have significant errors 12l-ll4l|. of order 40% or greater, 
we do not use them in exploring the constraints on X. 

An additional requirement which the X contributions to B K^*'>£~^£~ need to fulfill is that 
they do not upset the existing data on the processes B — )■ K^*'>cc followed by cc — £~^i~ in 
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TABLE I: Experimental branching-ratios of B ^ K'^*H+1- from Belle [13| and 

from CDF in units of 10~^, used to constrain the X contributions, for different ranges. The 

statistical and systematic errors have been combined in quadrature. 



;GeV2) 



[1,6] 
[14.18,16] 
> 16 



B{B Ki+i-) B{B+ K+fM+fi-) B{B KH+l') B{B^ ii:*V+/x- 



1.36 



i"0.24 
-0.22 



1.01 ±0.27 



1.49 
1.05 
2.04 



0.47 
0.42 
0.30 
0.27 
+0.31 
0.29 



1.60 ±0.56 
1.51 ±0.38 
1.35 ±0.39 



TABLE II: Standard-model predictions for branching-ratios of S — t- K^*H'^£ , in units of 10 ^, for 
different ranges, from Refs. j28|. 



(GeV2) 


B{B K£+£-) B{B KH+r) 


[1,6] 
[14.18,16] 
> 16 


J--Oo_Q 45 Z.DU_^ 34 

-I 00+0.43 
r-OZ_Q3g 

1 r;4+0.48 



the regions where the charmonium (cc) resonances lie, especially if m\ also falls within one 
of these regions. The relevant charmonia here are J/ip{l^) and ip' = V^(2S) whose masses are 
rrij^^ ~ 3.10 GeV and m^, ~ 3.69 GeV j22j, respectively. The constraints follow from comparing 
the experimental and SM values of the B — )■ {J/ip,ip')K^*\ {J/ip,ip') — )■ branching ratios, 
which are not available from the B — )■ K^*H~^£~ measurements, as events with bins in the 
neighborhood of Tn?]i^^j are vetoed to avoid these long-distance backgrounds mediated by the 



resonance (jj-uj. Experiments on B ^ J/^jK^*^ yield B{B+ J/^K+) = (10.14±0.34) x 10" 



B{B^ J/i)K^) = (8.71 ± 0.32) x 10-^ B{B+ J/i)K*+) = (14.3 ± 0.8) x 10-^ and -6(5° ^ 
J/i)K*^) = (13.3 ± 0.6) X 10"^ |22j, to be compared to the SM values B{B+ J/iIjK+) = 
(9.20;?;^^) X 10-^ i3(50 ^ J/^K^*) = (8.60lf:f ) x 10-^ BiB+ J/^K*+) = (9.95l^:?6) x 10"', 
and B{B^ J/^K*+) = (9.30^^|^) x 10"^ from Ref. (3. As for 5 ^ ^j'K, the data are 
B{B+ ^'K^) = (6.46 ± 0.33) x 10-^ i3(5° ^ ^'fsfO) = (6.2 ± 1.2) x 10-^ B{B+ ^'K*+) = 
(6.2 ±0.5) X 10~^ and B{B^ -> ■^'K*+) = (6.1 ±0.5) x 10"^ 0, whereas the only SM prediction 
of which we are aware is B{B^ — ip'K^) = 4.25 x 10~^ from Ref. [2^, without an error estimate. 
Since the SM calculations still have considerable uncertainties, after comparing the experimental 
and SM numbers for B{B -^{J/^,^')K^*^) B{{J/^,^') £+£-), with B{J/^ £+£-) ~ 5.9% 
and B{^jJ' — )■ £~^£~) ~ 0.77% [22], we may require 



-5 



-3x10-^ < B^{B^K£^£-)^,^^^^^^^^^^^^, < 5 x 10" , 

- 1 X 10-^ < B^{B^ ,..[8.6,10.2] GeV^ < ^ X 10"^ , (26) 

-1 X 10-^ < B^{B^ ^^*^^^^'),.e[12.8,14.2]GeV. < ^ >< 10"' ' 

where the chosen bins are similar to those for the charmonium backgrounds in the experiments. 
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To evaluate Bx{B — )■ K^*^t^l~')^ we employ the B — > K^*^ form factors from light-cone sum 
rules 26|, |27| . We have collected the formulas of the form factors in Appendix |Al In our search 
for the allowed parameter space of the X couplings, we use the lower bounds of the form factors 
in order to get the most space. 

After scanning the relevant parameters, we have found that there is X parameter space available 
which satisfies the constraints in Eqs. f l2T]) . fl2^ . ^IM . and fl2B]) . The size of the allowed parameter 
space can vary widely, depending on the mass and total width of X. To illustrate this, we take the 
X-mass values m-^ = 2, 3, 3.7, and 4 GeV. Since the couplings of X to other fermions are not 
specified in our approach, its total width is unknown. For definiteness we choose = 0.1 MeV 
in the four cases, but will comment later on other choices. Since in the decay amplitudes the 
couplings gy^ and g^^^ always occur each multiplied by gy^ or g^^, we show in Fig. [1] the allowed 
regions of their products in the cases of gy^ being real with g^^ = (top plots) and g^^^ being real 
with gy^ = (bottom plots). In the latter, gy^ = 0, cases, the restrictions from B — )• Ki'^i' are 
absent because g^^^ is not present in its amplitude. In the two rightmost plots, for = 4 GeV, 
the blue areas are invisible because they coincide exactly with the red ones. One could also get 
allowed parameter space for both gy^ being nonzero and one or both of them complex. We will 
have some examples with the complex couplings shortly. 

If one deals with gy^ separately from gy^^ instead of their products, it will be necessary to 
take into account constraints on the latter from the anomalous magnetic moments of the electron 



and muon, 



Og and a^. 



From their measurements, we have aT^^ 



115965218073 ± 28) x 10" 
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FIG. 1: Regions of (^y^, c/^^)Re c/^^ for Imgy^ 
9vs 



9 As = (top plots) and of {gy^,g^()Yliig^^ for 



Im^f^g = gy^ = (bottom plots) satisfying constraints from B — t- X^H.'^l^ (orange, lightly shaded), 
B Ki~^£~ (green, medium shaded), B — )• K*£'^i~ (blue, heavily shaded), and all of them (dark red). 
Prom left to right, the plots correspond to nix ~ ^' •^■'^^ ^ GeV. 
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and a?,^P 



exp _ 



(11659209 ± 6) X 10"^° 
= (-206 ± 770) X 10-^^ 



22| . The SM prediction for agrees with its measurement, 
but the SM prediction for presently differs by 3.2cr from its 



experimental value, a'j^^ — aT^ = (29 ± 9) x 10 [29|. Consequently, since the gy^^ and g^^^ terms 



in af are opposite in sign, we may impose 



- 9 X 10'^2 < < 5 X 10"^^ 



-1 X 10"^ < < 3 X 10" 



(27) 



Applying rux = 2, 3, 3.7, 4 GeV in the X contributions af in Eq. f|T7|) yields 



af(2GeV) = {5.5 g'y^ - 27. Q g^^) x 10-^° 
af (3GeV) = {2.Ag'y^ - 12.2 g^^) x IQ-'^ 



X 



(3.7 GeV) 



af (4 GeV) 



[1.6 g'y,- 8.1 gl,) x 10-^° 
(1.4(7^,-6.9(71) X 10-1° 



J(2GeV) = (22.8(7^^-117(7^J xlO-% 
(3 GeV) = (10.3 (7^^- 52.2 (7lJ x 10"^ , 
(3.7 GeV) = (6.8(7^^-34.3(7^^ x 10^^ , 
i;(4GeV) = (5.8 (7^^ -29.4 (7^ J x 10"^ . (28) 



X 



X 



To illustrate the parameter space of gyf> and (7^^ subject to the bounds in Eq. ( l27l) . we display in 
Fig. |2] the mx = 3 GeV case. One can conclude from Eq. ( l28l) and these plots that for each value 
of nix allowed {gv£:9Ai) region for £ = e is much larger than that for £ = /i, although the 



reverse is true for the imposed limits on af and in Eq. (|2 

We now explore the allowed gy^ ranges for gy^ values that can lead to predictions compat- 
ible, at the one-sigma level, with the DO anomalous measurement of the like-sign dimuon charge 
asymmetry in semileptonic 6-hadron decays [H, as we proposed in Ref. [2j. This implies that at 
least one of gy^ has to be complex. For rux = 2, 3, 3.7, and 4 GeV, using the results of Ref. jif, 
we obtain gy^ = (3.0 - 2.7i) x 10"^ (7.5 - 6.5i) x 10"^ (9.5 - 8.5i) x 10"^ and (12 - Hi) x 10"^ 
respectively, as possible values in the g^^ = case. In Fig. |3] we show the corresponding regions 
consistent with the various constraints described above. The areas displayed in the four plots have 
all turned out to be well within the bounds from both a, and (yellow, very lightly shaded, if 
not covered by the other bounds). In the fourth plot, for nix — 4 GeV, the blue area coincides 
exactly with, and hence is completely covered by, the red one. If gy^ = 0, one could also find the 
allowed parameter space, in which case the restrictions from B — )■ Ki^i~ (blue areas) are again 
absent and the regions consistent with all the other constraints are generally smaller than in the 
cases with g^^ = 0. 
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FIG. 2: Parameter space of {gve^dAe) ^^'^ {dvn^dAfi) '^x ~ 3 GeV subject to constraints from the 
anomalous magnetic moments of the electron and muon, respectively. 
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FIG. 3: Allowed ranges of gyg^ and for (from left to right) 



m 



X 



2, 3, 3.7, and 4 GeV with = 



and Qy^ values given in the text, subject to constraints from B — )• (orange, lightly shaded), 

B — )• Ki~^£~ (green, medium shaded), B K*£~^i~ (blue, heavily shaded), (yellow, very lightly 
shaded), and all of them (dark red). 



B. Predictions 



From the values of the X couplings gvs,As ^'^d gvi,Ae allowed by the various limits above, we can 
assess the effects that X might have on observables which will likely be studied experimentally in 
the near future. With improved precision compared to the existing data, the upcoming measure- 
ments will test the existence of X more stringently or place stricter constraints on its couplings. 
The observables we discuss here are the differential branching ratios of B — )■ K^*^i~^i~, the K* lon- 
gitudinal polarization fraction Fj^ in B ^ K*i~^i~, and the lepton forward-backward asymmetry 
Ap-Q in the latter decay, as well as the branching ratios of some rare decays. 

To illustrate how the X contributions may modify the SM predictions, we adopt for definiteness 
some of the larger values of the couplings from the top plots in Fig. [H Thus we have 



[dve^ 9ai)9vs 



m 



(1.0, -2.0) X 10"^^ for 

(-1.0, 0.5) X 10-9 for 
(1.0, -2.0) X 10-1° 

^ (-9.0, 3.0) X 10-11 for m 



= 2 GeV 



3 GeV 



for = 3.7 GeV 



4 GeV 



(29) 



with Qj^^ = 0. These numbers translate into the differential branching ratios of i? — ?■ K^*^i'^i^ 
in Fig. m where the curved (yellow) bands correspond to the SM expectation including the ±15% 
uncertainties in the form factors and for the curves corresponding to the combined SM and X 
contribution we have used the lower limits of the form factors. The two vertical (gray) bands mark 
ranges in which decay events are vetoed in the B — )■ K^*H^i~ experiments to reject backgrounds 



from the charmonium resonances 12-14 



From Fig. H] one can see that the impact of X can lead to mild changes in the SM branching 
ratios. In particular, the differential branching ratios may have spikes, indicating the X presence 



at q 



m 



Xj 



possibly accompanied by small dips and rises on opposite sides of the spikes arising 



from the enhanced interference near the resonant point between the SM and X amplitudes. To 
observe these features would require not only high precision in measuring the branching ratios in 
the bins, but also sufficiently small bin sizes. Both are not yet available in the existing data 



from BaBar, Belle, and CDF [iMA 
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FIG. 4: Differential branching ratios of S — t- Ki'^i~ (left plot) and B — t- K*i'^i^ (right plot) as functions 
of the squared dilepton-mass in the SM (yellow curved bands) and its combination with the X contribution 
for m-^ = 2 (green solid curves), 3 (blue dashed curves), 3.7 (red dot-dashed curves), and 4 (black dotted 
curves) GeV, with the gy^ and Qy^, values in Eq. 



With the couplings in Eq. fl29|) as before, we graph in Fig. |5] the K* longitudinal polarization 
fraction and lepton forward-backward asymmetry A-p-g for B — )■ K*i^i~. These plots suggest 
that, since the SM uncertainties are much reduced for these observables, the contributions of X to 
them, especially the latter, can produce modifications to the SM predictions that are likely to be 
more detectable than in the branching ratios. Thus the small spikes and deep narrow dips around 



9 



m 



corresponding to the dips and spikes in the differential branching ratios, could be very 
revealing. However, as in the branching-ratio case, it would be necessary to have high precision 
in the and A-p-Q measurements, as well as small sizes of the bins. This applies particularly 
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: Plots of K* longitudinal polarization fraction (left) and lepton forward-backward asymmetry 
for B — )• K*£^£~ in the SM (solid curves) and its combination with the X contribution for 

2 (green solid curves), 3 (blue dashed curves), 3.7 (red dot-dashed curves), and 4 (black dotted 
GeV, with the gy^ and Qy^ values in Eq. (f29]l . 
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to, for instance, the determination of the location of the Ap^ zero-crossing point which may only 
be slightly affected by the presence of X, as the right plot in Fig. [5] indicates. In contrast, the 
available experimental results on and A-p^ still have substantial relative errors, more so than 
the branching ratios 



12 



-14 



At this point we should comment that the effects of X as depicted in Figs. S] and correspond 
to the choice for its total width Fx = 100 keV which we made earlier. If a larger (smaller) value of 
were assumed, the parameter space of X consistent with the constraints we discussed would be 
larger (smaller) also, and the signals of X would become more (less) detectable. In any case, the 
availability of the allowed parameter space should encourage experimental searches for X in future 
6-hadron experiments. 

Finally, we turn to a few rare decays, the first one being — )■ , which also proceeds 

from b — )■ si~^£~. The expressions for its amplitude and decay rate can then be obtained from 
those for B — > K*i'^i~, with appropriate changes. Using the -B^ — )■ form-factors listed in 
AppendixEl we estimate the SM branching ratio to be 1.0x10"^ < Bg^{B^ (f)£+£~) < 1.8x10"^ 
after excluding 8.68 GeV^ < < 10.09 GeV^ and 12.86 GeV^ < < 14.18 GeV^ from the 
integration, as in the experimental studies. Very recently CDF has reported the first observation 
of this decay, with B^^p{B^ (pf^^f^~) = (1-44 ± 0.33 ± 0.46) x lO'^^ [U], which is compatible 
with the SM value. In view of the still-large theoretical and experimental errors, there is room 
in this decay for the X contribution, which can expectedly be tested with improved precision 
at LHCb. We further expect that the longitudinal polarization fraction and lepton forward- 
backward asymmetry in this decay are sensitive tools to detect the signals of X, much like their 
counterparts in i? — )■ K*£~^i~ . 

As mentioned earlier, the decays B^ — > {ri,ri')i~^£~ are also relevant in the search for X, but 
they are yet to be observed. They may be detected in the near future, as the SM values of their 



branching ratios are estimated to be only several times smaller than those of B^ — )• (f)i~^i~ [31 
Since the amplitudes for B^ — )■ {ri,rj')i'^i~ are independent of the coupling g^^, as in the case of 
B — ?► K£~^£~ , we expect that qualitatively the effects of X on the former would be roughly similar 
to that on the latter. 

Additional decays of interest which are not yet observed are B^ — )■ Their amplitudes 

do not involve gy^, in contrast to B^ {ri,ri')£~^i~ . The experimental upper-bounds on their 



branching ratios are available 22 



B,^^{B, ^ e+e-) < 2.8 x 10^^ , B,^^{B^ ^ < 3.2 x 10"^ . (30) 

From Eq. (fT6l) . we find the SM values 



BsM{Bs^e+e-) ^ 7.5xl0~'\ Bsu{B, ^ fi') ~ 3.2 x 10"^ (31) 



using the central value of Jb^ = 240±30MeV 3(|. If we include the contribution of X in Eq. fll6p 



we have, for instance, —9 x 10 ^ dAsdAi ~ 6 x 10 from the second of the bottom plots 
in Fig. [H which translates into 



-14 



6.6 X 10"^^ < B{B^ e+e") < 8.1 x 10" 
2.8 X IQ-^ < B{B^ fi+ii-) < 3.5 x 10"^ . (32) 
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Hence the X contributions arc easily accommodated by the present experimental limits and can 
produce moderate modifications to the SM predictions. Once Bg — )■ fi'^fi~ has been seen, it 
can provide important restrictions on the combination QasOah branching ratio turns out be 

consistent with the SM expectation. 



IV. CONCLUSIONS 



We have extended our proposal in an earlier paper that a nonstandard spin-1 particle lighter 
than the h quark with flavor-changing couplings to h and s quarks can offer a viable explanation for 
the recent anomalous measurement by DO of the like-sign dimuon charge asymmetry in semileptonic 
6-hadron decays. Specifically we have considered the possibility that the new particle also couples 
to the light charged leptons £ — e, fj, and thus contributes to rare b ^ s transitions involving 
the leptons. After exploring experimental constraints on its couplings from currently available 
experimental data on the inclusive B Xsi^i~ and exclusive B K^*^i~^i~ decays, as well as 
the charmonium contributions to the latter and the anomalous magnetic moments of the leptons, we 
have shown that there is parameter space of the particle that is consistent with the data, including 
the DO result. This can serve to help motivate dedicated searches for X in &-hadron experiments in 
the future. We have found that the differential branching ratios of these decays, the K* longitudinal 
polarization in B — )■ K*i^i^, and its lepton forward-backward asymmetry, especially the latter 
two, are observables that are potentially sensitive to the signals of X. Furthermore, rare decays of 
the Bg meson, such as Bg {((),r],7]')i'^i~ and Bg — >■ can provide additional observables. 

We expect that the upcoming measurements of rare B and Bg processes with greatly increased 
precision at LHCb and next-generation B factories can probe the existence of the particle, or its 
couplings, stringently. 
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Appendix A: Form factors 

To obtain the B — )■ Ki'^£~ amplitude, we use the matrix elements 
{K\sYb\B) = 

{K\sb\B) = !<^F, 



— 



) 



/-I I -\ {Pr + PkY — i^'n — TTl^c^Q^ 
qMsaf'^B) = ^ ^ir^iF^^ 

{K\s^^^,h\B) = {K\s^,h\B) = q''{K\-sa^^^,h\B) = , (Al) 
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where {pj^) is the B (^K) momentum, q = — Vk^ ^^"^ form factors Fg ^ each depend 
on q^. The matrix elements relevant to i? — t- K*i 



, are 



[K*\-siAb) 



2V 



e*-q 



q 

iA^ e*-q 
nig + iT^K* 



q^ 



[K*\sh\B) = 0, {K*{pj,,)\s^,b\B{pg)) 



-2iAQ rrij^* e*-q 



(f{K*\s(r^^-i^h\B) = T2 {ml-ml,)^^- {pg+pj^.) e*-q 



iPB+PK*)f,Q^ 



e-q 



(A2) 



where pp^, is the K* momentum, q = p^ ~ Pk*i ^^"^ form factors V , A^ ^ 2-1 ^^^d Tj^ 2 3 
functions of g^. 

In numerical calculations, we adopt the form-factor results of Refs. 26|, \2j\ from light-cone sum 
rules. The form factors are parametrized as 



F{s) 
F(s) 



for F^, 



F(s) = %^ for Fq , , T2 , 



1 - s/ml^ 



+ 



l-s/m\ 1-sMfit 



+ 



l-^MIt [l-s/ml,]' 



for V, A,, T, 



for A^, T3 , 



(A3) 
(A4) 
(A5) 
(A6) 



where T3 = + sT^^jim^ — my), with = m^o.mQ^ and my = mj^*,m^. The parameter 
values in the B — )■ R'^*^ and B^ ^ (f) cases are collected in Tables [TTll and HV] respectively. The 
estimated uncertainty of each F(s) is about ±15% |26|-|28 



TABLE III: Parameters for B — )• K^*^ form factors 26|, |27|. The and entries are in units of GeV'^ 
and GeV, respectively. 









Ft 


^0 






V 


Ti 


T2 


^3 


^1 




0.162 


0.161 


1.364 




-0.084 


0.923 


0.823 




-0.036 


^2 


0.330 


0.173 


0.198 


-0.990 


0.292 


0.342 


-0.511 


-0.491 


0.333 


0.368 




37.46 






36.78 


40.38 


52.00 


49.40 


46.31 


41.41 


48.10 






5.41 


5.41 


5.28 






5.32 


5.32 
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TABLE IV: Parameters for ^ <j) form factors [27]. The m|^. and entries are in units of GeV^ and 
GeV, respectively. 







^1 


^2 


V 




T2 


^"3 


^1 


3.310 




-0.054 


1.484 


1.303 




0.027 


^2 


-2.835 


0.308 


0.288 


-1.049 


-0.954 


0.349 


0.321 




31.57 


36.54 


48.94 


39.52 


46.31 


37.21 


45.56 




5.37 






5.41 


5.41 







Appendix B: Squared amplitudes and decay rates 

The general expressions for the (double) differential decay rates of -B — )■ K^*H~^i~ arising 
from various b — )■ si~^i~ four-fermion operators within and beyond the SM are known in the 
literature [l^, [13] ■ Here, for completeness, we write down the specific formulas resulting from our 
amplitudes of interest. 

The absolute square of the B — )■ K£'^i~ amplitude in Eq. (Q, summed over the lepton spins, is 
\M{B ^ Ke+e-)\' = ^^^M^I\A\'[{t-m^){u-m^)-mlml] 



+ |Cp [tu~ {ml - m^) {mj, - mj)] 
+ \D\^mlq^ + 2Re{C*D) ( 



m R — ) m 



(Bl) 



with t = {j>^ — Pi+Y ^ ~ {Pb ~Pe-Y- '^^^^ leads to the differential decay rate 



210 



^^i\\A?[i-'-^\^\C\ 



e-^ + (2 + 2p-z)(l 



+ |Dp(l-i;2)z2 + 2Re(C*D)(l-p)(l-t;2);2} , (B2) 



where ^ = 1 — 2p — 22; + (p — 2;)^, p = m^/m^, v = a/1 — Amj/q'^, and z = /m\. 

The absolute square of the B — )■ K*i~^£~ amplitude in Eq. ( !TT|) . summed over the K* polariza- 
tion and lepton spins, can be written as 

„,2r'2 I \ 12^2 
\M{B ^ K*i^r)\' = " X 



24 TT^r 



+ |£p)m| [a(2 - v') + [i - uYj r z + 4|C|V(3 -v^)z + 8|^ 



1 2 2 
rv z 



+ 



\C\^ + \J^?+ {\V\^ + \Q\^)\m% + 2Re{C*V + rg)ml{r + z-l) A- (t - m) 
2|£|V + \g\^{z - 2 - 2r) - \H\'^z + 2Re{g*H) (r - 1)1 Am^(l - v^)z 
- 2 Re{J^*g + J^*H) A m|(l - 2 + 8 Re(^*J^ + C*£) r (t - m)^} , (B3) 
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where \ = 1 — 2r — 2z + (r — z)"^, r = m|-./m^, t = t/ni^, and u = u/m\. The corresponding 
differential decay rate is 



212 



2\AYrz + ^ 1 + 



B 



3 r 
12r;z 



X 



A 



(3-t;2) + 4|£|Vt;2z + 3 |7/|2(i _ 



+ 



24r 2 



A 



+ \g\^ (1 - r)2(3 - v^) + 2(z - 2 - 2r) 



+ 



2Re(C*P), 2\ 2Re(J'*^) 



+ 2- l)(3-t;2) + 



B 



(1 -r)(f2-3) +2i;2z 



+ 6 



Re(J^*'H) 



Re(e;*H) (r - 1) 



B 



(B4) 



From Eqs. ( fT3l) and ( ]B3[) . we can derive the i^'* longitudinal-polarization fraction 



ml 



+ Re(e;*-H) (1 - r) 



+ 



rrX 



B 



, nx z 1 + 
2r(3-t;^)- + ^- 



+ 



+ 



PpA Re(C*P) 



+ 



(r + 2 - 1) 



(1 + .^) 



+ — ^ 



A 



+ (1 - r)' 



+ ^^^?^[r-l + (r + 2^-l)i;2] 



(B5) 



and lepton forward-backward asymmetry 

1 



A 



FB 



1 d^r(B ^ K*e+r) 

dca ^ — r^r^ sgn 



(ig2 (icg 
Re(^*J^ + C*^) . 



(B6) 



2107^5^/(5 ^ K*i+i-) 

where r'(5 KH+t) =dT{B ^ K*£+£-)/dq^ and = cos0. 

We remark that the t — u term in the last line of Eq. (IB 3 1) contains Apg for B — >■ K*i^i~ . 
Since such a term is absent in Eq. (IBip . there is no Ap^ for 5 — > Ki^i~ from the SM or from the 
X contributions under consideration. 
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